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Abstract 

We study nondeterministic and probabilistic versions of a discrete dy- 
namical system (due to T. Antal, P. L. Krapivsky, and S. Redner [3J) in- 
spired by Heider's social balance theory. We investigate the convergence 
time of this dynamics on several classes of graphs. Our contributions 
include: 

1. We point out the connection between the triad dynamics and a gen- 
eralization of anmhilating walks to hypergraphs. In particular, this 
connection allows us to completely characterize the recurrent states 
in graphs where each edge belongs to at most two triangles. 

2. We also solve the case of hypergraphs that do not contain edges 
consisting of one or two vertices. 

3. We show that on the so-called "triadic cycle" graph, the convergence 
time is linear. 

4. We obtain a cubic upper bound on the convergence time on 2-regular 
triadic simplexes G. This bound can be further improved to a quan- 
tity that depends on the Cheeger constant of G. In particular this 
provides some rigorous counterparts to experimental observations in 

m- 

We also point out an application to the analysis of the random walk 
algorithm on certain instances of the 3-XOR-SAT problem. 

Keywords: social balance, discrete dynamical systems, rapidly mixing 
Markov chains, XOR-SAT. 

1 Introduction 

Discrete Dynamical Systems, both synchronous (a.k.a. cellular automata [19]) 
and asynchronous (e.g. sequential dynamical systems [22]) provide a rich family 
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of models, whose computational properties (e.g. computational universality, the 
complexity of prediction, or properties of the phase space, such as the existence 
of Garden of Eden states) have been actively investigated. Such models under- 
line some of the most celebrated examples in the social sciences (e.g. the cel- 
ebrated Schelling Segregation Model) and have, unsurprisingly, been proposed 
(see e.g. [TB], U) as the basis for modeling social dynamics. 

In this paper we are concerned with the properties of a particular dynamics, 
due to T. Antal, P. L. Krapivsky and S. Redner [3]. The dynamics (referred in 
the sequel as the triadic dynamics) is inspired by Heider's balance theory [17], a 
well-established subject in social psychology, and one of the first theories in the 
social sciences to benefit from concepts and methods from graph theory [5]. 

Heider's balance theory aimed to describe the "equilibrium" properties of the 
interpersonal relations, and did not include a dynamical component. Recently, 
however, dynamical models related to the one we study, have become popular 
in the statistical physics [211 EO] and social simulation literature [HI |29] . 

The triadic dynamics (described precisely in Section 2 below) is parametrized 
by a constant p € [0, 1]. In [^ the authors investigated the dynamics on the 
complete graph Kn and displayed a double phase transition in the convergence 
time with respect to control parameter p. The triadic dynamics was further 
investigated by means of computer simulations in [25j . The graph topology in 
this paper is a finite section of the triangular lattice. A phase transition is also 
displayed around a critical value pc ~ 0.4625. In this case the convergence time 
is polynomial for p < pc and logarithmic for p > pc- 

We will be concerned with extending the study of the triadic dynamics to 
general graphs. To do so, we have to solve a problem that has appeared many 
times in the study of dynamical systems, that of state reachability: when is a 
given system configuration S2 reachable from a given initial state si by a finite 
sequence of moves ? In general (e.g. [6l|5]) this type of problem is computation- 
ally intractable even for very simple dynamical systems. In contrast our main 
result shows that for a large class of graph topologies questions such as reach- 
ability are computationally tractable for our dynamics. We complement this 
result with a study of the convergence time of our dynamics in the probabilistic 
setting similar to that of [3]. 

A completely different reason for our interest in the triadic dynamics is 
its somewhat unexpected relation [26] to 3-XOR-SAT, a combinatorial prob- 
lem that was investigated in Theoretical Computer Science [11] and Statistical 
Physics [27]. In particular, using the methods we develop for the analysis of the 
dynamics we will upper bound the expected convergence time of a simple local 
search algorithm, called Random Walk, previously analyzed for random 3-XOR- 
SAT in [28l|7]. Unlike these papers, we will analyze the convergence time of the 
algorithm on individual formulas, and bound the convergence time (Theorem [5] 
below) using quantities that depend on the structure of the input formula. 
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2 Preliminaries 

The following is a formal definition of the dynamics we will be investigating in 
this paper: 

Definition 1 Nondetermininstic Triadic Dynamics. We start with a graph 
G = {y,E) whose edges are labeled ±1. A triangle T is G is called balanced 
if the product of the labels of its edges is equal to 1. At any step t, for any 
imbalanced triangle T we are allowed to change the sign of an arbitrary edge of 
T (thus making T balanced). The move might, however, make other triangles 
unbalanced. 

Definition 2 Probabilistic Triadic Dynamics. A probabilistic version of 
the dynamics in Definition]^ parametrized by a real number p £ (0, 1), is spec- 
ified l&j as follows: While there exists an imbalanced triangle, first choose uni- 
formly at random an imbalanced triangle T . If the triangle T has a single 
negative edge e: 

• With probability p turn e to positive. 

• With probability 1 — p turn one of the other two edges of T to negative. 

otherwise (T has three imbalanced edges ) change the label of a random edge of 
T. 

Definition 3 A triadic simplicial complex is a graph G — (V, E) such that all 
edges e E are part of some triangle of G. 

Definition 4 The triadic dual of graph G is an undirected hypergraph with 
self-loops Tz{G) = (1^3, i^s) defined as follows: is the set of triangles of G. 
Hyperedges in T'i{G) correspond to edges in G and connect all vertices v e V3 
containing a given edge. In particular we add a self-loop to vertex v € V3 if v 
belongs to an unique triangle of G. We may even add two self-loops to the same 
vertex T if two of its edges belong only to triangle T . 

We will also require that the triadic dual of graph G be connected. If this 
was not the case, the dynamics would decompose on independent dynamics on 
the connected components. 

Definition 5 The triadic cycle TC„ is the graph G consisting of n triangles 
chained together. Edges AB and CD of the extreme triangles in the chain are 
"glued" by identifying A = C , B = D. 

The construction is illustrated in Figure 1, where the triadic cycle graph 
with sixteen triangles is displayed, together with its triadic dual (in this graph 
circles represent self-loops). 
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Figure 1: (a). The triadic cycle TCis- (b). Its triadic dual 



Definition 6 Let G = {V,E) he a graph. For S C V denote by vol{S) = 
SxGS d{x) and E{S, S) the set of edges connecting a node in S to one in S . 
The Cheeger time of graph G is 

Tc[G) = sup =^ =^ f- , 

where 'n[A\ = ^ and, fori,j G V{G), we define pi j — l/deg(i) iff i and j 
are adjacent, otherwise. 



Remark 1 WhenG is r -regular Tc{G) = sup^^^ 



r/n-\E(A,A)\ 

Finally, we need the following concept from [T]: 



sup^cv 



f rn-\E(A,A)\ 



Definition 7 Let (F„) be independent random walks on G The worst case 

meeting time of G is defined as: 

tm{G) = max E[Tm\Xo = = j], 
«,iev(G) 



where Tm = min{fc : Xj^ = Yfc}. 



3 Connection with annihilating walks 

The triad dynamics can be mapped to a generalization of annihilating walks. 
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Definition 8 An annihilating walk (AW) on graph G is a nondeterministic 
dynamical system described as follows: Start with a ball at some of the vertices 
of the random graph. At each step we move one ball to a neighbor of its current 
vertex. When two balls meet at a vertex v they annihilate each other, and are 
subsequently eliminated from graph G. 

The stochastic version of annihilating walks (called annihilating random 
walks [13]) have been studied in the interacting particle systems literature ( 
e.g. [M]). A few recent results [10], [2] deal with interacting particle systems on 
a finite graph as well. 

We will also need the following: 

Definition 9 A coalescing random walk (CRW) on graph G is a stochastic 
process described as follows: 

1. Start with a ball at some of the vertices of the random graph. 

2. Each ball performs a random walk on G. 

3. When two or more balls meet at a vertex v they coalesce, performing a 
common random walk on G. 

Definition 10 For an annihilating random walk we denote by TABwiG) the 
random variable defined as the number of steps until the number of balls on graph 
G becomes zero (one if the initial number of balls was odd). For a coalescing 
random walk we denote by Tcrw{G) the random variable defined as the number 
of steps until all the balls coalesce. 

To recover the connection between annihilating walks and triad dynamics, 
we have to define a generalization of annihilating walks to hypergraphs. 

Definition 11 The hyperedge switching process (HS) on a hypergraph G is 
a nondeterministic dynamical system described as follows: Start with a ball at 
some of the vertices of the hypergraph. At each step we choose a hyperedge C 
in G that contains at least a ball. Any self-loop on a vertex that contains a ball 
counts as one of these such cliques. We then put balls on all empty nodes of C 
and remove the balls from the nonempty nodes. 

Proposition 1 The nondeterministic triad dynamics on graph G and initial 
configuration s corresponds to the hyperedge switching process on graph T^{G) 
to the initial configuration s defined by putting a ball to all vertices v £ T3{G) 
corresponding to unbalanced triangles in G. 

Proof. This is immediate by duality. Indeed, choosing an edge e to flip (in the 
original graph) corresponds to chosing the corresponding hyperedge le (in its 
dual hypergraph). The choice is only possible if there exists at least one imbal- 
anced triangle containing e, that is (by duality) at least one of the vertices of 1^ 
contains a particle. Flipping the sign of edge e only affects triangles containing 
it (i.e. vertices of le in the triadic dual), and has the effect of making imbal- 
anced triangles balanced and viceversa. This corresponds plainly to removing 
balls from the occupied vertices of le and adding balls to the free vertices. (} 
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4 Recurrent states for nondeterministic triad dy- 
namics 

Definition 12 A state s is specified by giving a ±1 label to each edge of graph 
G. A state is recurrent if for any infinite path tt = (ttq, . . . , 7r„ . . .) of the system 
s is reachable from iTn ■ 

It is fairly easy to determine recurrent states if every edge is part of at most 
two triangles (that is, if the triadic dual is a graph with loops): 

Theorem 1 Let G be a graph such that (i) the triangle graph of G is connected, 
and (ii) every edge of G is part of at most two triangles, and let sq : E{G) — > 
{±1} be an initial state. Then: 

1. If G contains an edge that is part of only one triangle, then the recurrent 
states for the triad dynamics are exactly the completely balanced states 
reachable from sq. 

2. If every edge in G is part of exactly two triangles and the number of 
imbalanced triangles in Sq is even then the recurrent states for the triad 
dynamics are exactly the completely balanced states reachable from sq ■ 

3. If every edge in G is part of exactly two triangles and the number of 
imbalanced triangles in sq is odd, then the recurrent states for the triad 
dynamics are the states containing exactly one imbalanced triangle. 

Proof. 

Define function P that maps each state s to the number of imbalanced 
triangles in state s. 

Proposition 2 Function P is nondecreasing. That is, if s ^ t is a legal tran- 
sition of the dynamics, then P{s) > P{t). 

Proof. Consider the edge e G G chosen by the dynamics. There are three 
cases: 

1. There is only one (imbalanced) triangle T that contains e. Then, by 
flipping the sign of edge e triangle T becomes balanced and no other 
triangle is affected. The number of imbalanced triangles goes down by 
one. 

2. There are two triangles Ti,T2 containing edge e (Figure [3]), both im- 
balanced. Then, by flipping the sign of edge e triangles Ti,T2 become 
balanced and the balancedness of no other triangle is affected. Hence the 
number of imbalanced triangles goes down by two. 

3. There are two triangles Ti, T2 containing edge e, one imbalanced (say it is 
Ti) and one {T2) balanced. Then, by flipping the sign of edge e triangle 
Ti becomes balanced, while T2 becomes imbalanced. The balancedness 
of no other triangle is affected. Thus the number of unbalanced triangles 
stays the same. 
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Consider now a recurrent state t. It follows that there is no state w reachable 
from t with P{w) < P{t). 

To prove 1 note that from any state containing some imbalanced triangle 
we can reach a state w with one less balanced triangle by "propagating the 
imbalance" towards the triangle T with an edge e occurring only in T. In this 
process the number of imbalanced triangles stays the same. By then changing 
the label of e we decrease the number of imbalanced triangles. 

The arguments for points 2 and 3 are similar, noting that the number of 
imbalanced triangles stays the same or decreases by exactly 2. () 

We would like to extend our result to general graphs G. However, this is an 
open problem so far. The reason is that for general graphs G the triadic dual 
is no longer a graph, but a hypergraph and the potential function is no longer 
nonincreasing. In the sequel we provide a partial result. 

4.1 Recurrent configurations for the nondeterministic hy- 
peredge switching process 

In this section we study the hyperedge switching process on hypergraphs H 
whose hyperedges all contain more than two vertices. We will refer to this 
condition as H does not contain graph edges. 

Definition 13 For every pair of boolean configurations wi,W2 ■ V{H) — > Z2 
on hypergraph H we define a system of boolean linear equations H{wi,W2) as 
follows: Define, for each hyperedge e a variable Zg with values in Z2. For any 
vertex v G V{H) we define the equation 

= W2(w) - (1) 

In equationUlthe difference on the right-hand side is taken in Z2; also, we allow 
empty sums on the left side. System II{wi,W2) simply consists of all equation 
Qj, for all V e V{H). 

Definition 14 If x is a state on H and I is an edge of H , define 

I x(v), otherwise 

In this section we prove the following: 

Theorem 2 Let H be a connected hypergraph with no graph edges, let wi be an 
initial configuration that is not identical to the "all zeros" configuration 0, and 
let W2 be a final configuration. 

There is a polynomial time algorithm to test whether W2 is a recurrent state 
for the nondeterministic hyperedge process on H with starting state wi . 

Proof. 



(2) 
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Proposition 3 // state W2 is reachable from wi then the system of equations 
H{'Wi,W2) has a solution in Z2. 

Proof. Let P be a path from wi to W2 and let Ze be the number of times edge 
e is used on path P (mod 2). Then {ze)eGE is a sohition of system Indeed, 
element w{v) (viewed modulo 2) flips its value anytime an edge containing v is 
scheduled. (} 

Lemma [3] gives a necessary condition for reachability in Theorem [2] The 
complete characterization of recurrent states is a consequence of the following 
two lemmas: 

Proposition 4 Let H be a hypergraph with no graph edges (not necessarily 
connected). Assume that wi,W2 G {Ojl}^*-^^ are configurations such that for 
no connected component C of H, wi|c = 0, W2\c ^ 0. Then state W2 is 
reachable from wi if and only if system of equations H{wi, W2) has a solution 
in Z2. 

Proof. We prove the result by induction on the number of edges in hypergraph 
H. 

1. Case m = 1: Suppose system H{wi,W2) has a solution. Since H con- 
tains a single edge e, W2(v) = wi{v) for all vertices v ^ e (otherwise the 
system would contain equation = 1). Also, for all w G e the quan- 
tity W2{v) — wi{v) := A does not depend on v (otherwise we would have 
two contradicting equations). There are two cases: A = 0, in which case 
wi — W2 and path P can be taken as the empty path, and the case A = 1. 
In this case it follows that edge e contains at least one vertex v with 
'Wi{v) = 1. Indeed, if this was not true then ifije = 0, and W2\e = 1, con- 
tradicting the restriction from the hypothesis. One can then set P = {e} 
and obtain a path from wi to W2. 

2. Case m > 2: 

Assume Lemma |4] is true for all hypergraphs with less than m edges, 
and consider an arbitrary hypergraph H with m edges. Without loss of 
generality we may assume that H is connected, otherwise we find paths 
separately on every connected component. There are two cases: 

(a) For some I G E system H{wi,W2) has a solution with z; = 0. 

Then the system U{wi, W2), corresponding to hypergraph U — H\{1} 
is solvable; indeed every solution of H{wi,W2) with zi — is also a 
solution of U{wi, W2). 

Hypergraph U has four types of connected components: 

i. Connected components P oiU such that t«i|p\{;} ^ 0. 

ii. Connected components Q of U such that wi\q = but there 
exists V €z I r\ Q with wi(v) — 1. 

iii. Connected components R of U such that wilu = 0, W2\r =/= 0. 
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iv. Connected components S U such that wi\s = 0, 'W2\s = 0. 

Since H was connected, all such components contain at least one ver- 
tex from I. Moreover, we can assume that there exist no components 
of type (iv) , since we can eliminate them from consideration as they 
do not affect the overall result. Also, the induction step is trivial if 
U has no components of type (iii). So assume that U contains some 
component of type (iii). 

• Case 1: wi\i ^0 and there exists z £ / such that either 
Wi(z) = cind z belongs to a component of type (i) or 
W2{z) = and z belongs to a component of type (iii). 
We update the state using the following sequence of steps. 

— 1: Schedule edge I. This is possible since the current con- 
figuration contains, by the hypothesis, at least a node v € I 
with wi{v) = 1. In turn, this changes the states of nodes in 

I that belong to components of type (iii) to one. 

— 2: Schedule components of type (iii), moving the states of 
their nodes from to ^2^. This is possible since systems 
U{wi,W2) and U{'Wi\'W2^) are equivalent (thus satisfiable), 
and state contains a nonzero value on every connected 
component of type (iii). 

— 3: Schedule edge I again. This is possible since the state 
of node 2: is 1. This makes the states of nodes t in I that 
belong to components of type (iii) take value W2{t). Also, 
this changes the values of other nodes t €l to wi [t] again. 

— 4: Finally, schedule components of type (i) and (ii), changing 
the values of their nodes from wi to W2- This is done com- 
ponentwise, using the induction hypothesis, and the nonzero 
values of wi in each such component. 

• Case 2: wi\i ^ 0, Case 1 doesn't apply and there exists 
P a component of type (i) or (ii) with ^2(2) = 1 for some 

zGlnP. 

Since Case 1 docs not apply, we can identify the values of wi 
on every vertex x & I: wi{x) = 1 if a; belongs to a connected 
component of type (i) or (ii), wi(a;) = otherwise. Also, if x G I 
belongs to a configuration of type (iii), W2{x) = 1. 
Consider an edge I2 such that I2 r\ I 7^ and I2 belongs to a 
component of type (iii). I2 exists (i.e. Q does not consists of an 
isolated vertex in U) otherwise the component Q would yield an 
unsatisfiable equation = 1 in U{wi,W2)- 

— 1: Schedule components of type (i) and (ii), changing the 
values of their nodes from wi to W2- 

— 2: Schedule edge I. This is possible since the current con- 
figuration contains, by the hypothesis, at least a node z G I 
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with 'W2{z) = 1. In turn, this turns the states of nodes in / 
that belong to components of type (iii) to one. 

— 3: Schedule all components P of type (iii) except Q, mov- 
ing the states of their nodes from wj'-* to 'W2 \ restricted to 
the scheduled components. This is possible since systems 
U{wi,W2) and U{wi \w2 ^) are equivalent. 

— 4: Schedule component Q, moving its state from w^^Iq to 
{w2^)\q. This is possible since system U{wi\w2^) is solv- 
able. 

Note that, since Case 1 doesn't apply, W2 \x) = for all 

X eQnl. 

— 5: Let r be a node of Q with 'W2{r) = 1, closest to some 
vertex of edge I in distance. Consider the path D from r to a 
node in I. All edges other than the one containing r have only 
zero values. Schedule edges on path D from r towards I, one 
by one (Figure 3(a)). This "propagates" the value 1 towards 
the vertex in I. Note that all nodes on intermediate edges 
had initially zero values, so they have intermediate nodes on 
the path retail value 1. In particular, last edge has a node in 
Q\{1} whose value is one. 

Also, this action turns the values of nodes in ^ n Q to one. 

— 6: Schedule edge I again. This restores the correct value of 
all nodes in l\Q. It also turns nodes in l\Q to zero. 

— 7: "Undo" scheduling nodes on D from / to r (Figure 3(b)). 
This is possible since, as each edge has > 3 vertices, last 
edge had a node whose value is one, and "we propagate this 
value" towards r. Since edges of D were scheduled twice, 
every label is correct. 

• Case 3: wi\i ^0 and Cases 1 and 2 do not apply. 

Since Cases 1 and 2 do not apply, we can identify the values 
of Wi,W2 on every vertex x € I: wi{x) — l,W2{x) = if a; 
belongs to a connected component of type (i) or (ii), Wi{x) = 
otherwise. Also, if x € I belongs to a configuration of type (iii), 

'Wi{x) = 0, «;2 (•'/■) = 1. 

The construction is almost identical to that of Case 3. The only 
difference is that Step 1 is now executed at the very end of the 
process. This is possible, since before executing this step all 
vertices x € I belonging to components of type (i) or (ii) have 
value 1 (edge I being scheduled twice was the only action affecting 
their value). 

• Case 4: wi\i = 0. We reduce this case to one of the previous 
three cases as follows: since wi is not equal to 0, there exists a 
vertex v reachable from / in H with wi{v) = 1. Choose such a 
vertex v at minimal distance from I . Let this path be D. First, 
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Figure 2: (a). Forward propagation of ones. (b). Backward propagation. In 
both cases edge I is shaded. 

use the forward propagation trick on D to change the state wi to 

Cl StrCltrG u>3 with w^li 7^ 0. System H{w3,'W2) = has a sohition 
with zi = 0, since system H(wi,W2) = has one. We then apply 
one of cases one to three, the one that works, to the former 
system. 

(b) For all edges I, all solutions of system H{wi,W2) = have 
zi = 1. 

If there exists I with 'W\\i ^0 and twf ^ ^ then system H{wi^ ,W2) = 
has a solution with zi = (any solution of H(wi, W2) = with the 
value of zi flipped to zero). Wc then apply one of the previous cases. 

Otherwise for all I with wi\i ^ 0, Wi '' = 0. Therefore H has a single 
such edge, and the result immediately follows. 

Proposition 5 Let H be a connected uniform hypergraph, and let ■wi,'W2, W3 be 
configurations on H. If W2 is reachable from wi and ws ^ is reachable from 
Wi in one step then W2 is reachable from w^. 

Consequently one of the following two alternatives hold: 
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• is reachable from wi and is the only recurrent state for the process started 



• is not reachable from Wi and W2 is a recurrent state for the nondeter- 
ministic hyperedge switching process on H started at wi . 

Proof. If system H{wi,W2) has a solution (z^) and I is the edge flipped when 
going from wi to wg, let 



It is easy to see that (z^) is a solution to system H{w3, W2)- Since W2 is special, 
the restriction in Lemma [4] is satisfied for the pair {w^, W2)- <C> 

Applying Lemma|4]and Lemma[5]we obtain an easy algorithm for recurrence: 
make two reachability tests using Lemma ID <^ 

5 Time to social balance on the triadic cycle 

Theorem [T] motivates the study of the following quantity: 

Definition 15 For a triadic simplicial complex G, denote by tsb{G) the max- 
imum, over all initial states s, of E\Ts(G)], where Ts{G) is the expected time 
for the system, starting from initial state s to enter a state with zero (one) im- 
balanced triangles (according to the result in Theorem]^. tsb{G) will be called 
the time to social balance of graph G. 

We first investigate the convergence time of random triad dynamics on the 
triadic cycle T!„. For this class of graphs, Theorem [1] states that all recurrent 
states are socially balanced. We were motivated in our particular choice of graph 
G by our paper [12j . which investigated a similar dynamics. The dynamics in 
[T^ can also be described using the particle analogy, and is specified by particle 
creation/annihilation rules (1,0) -> (1,1), (0,1) (1,1), (1,1) (0,0). In 
contrast, the annihilating random walk corresponding to the dynamics in this 
paper is specified by a very similar set of rules (0,1) — » (1,0), (1,0) (0,1), 



The next result shows that the convergence time to a socially balanced state 
is similar to the analog result in [12]. The proof is simpler, though: 

Theorem 3 Let p^ — min{p, 1/2}. The time to social balance on the triadic 
cycle Tn isO{ f.). 

Proof. 

The result follows from the following: 



at wi . 




Ze, if e 7^ 

1 + z;, otherwise. 



(3) 



(1,1)^(0,0). 



Proposition 6 Function P decreases at each step with probability at least p_ 
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Proof. 

Consider three adjacent triangles in the graph TC„, as displayed in Figure |4l 
Suppose that triangle B is unbalanced and is the one chosen by the dynamics. 
There are two cases: 

1. The dynamics changes the value of edge x (the one that is not part of any 
other triangle). In this case no triangle but B is affected, and the number 
of unbalanced triangles goes down by one. 

2. The dynamics changes the value of one of the edges y, z. In this case 
triangle B becomes balanced. The triangle that shares the chosen edge 
with B becomes unbalanced if it was previously balanced, and balanced if 
it was unbalanced as well. All in all, the number of unbalanced triangles 
does not increase. 




Figure 4: The three triangles in the 
proof of Lemma [6l 





We now apply the following result ([24] Theorem 1.3 pp. 15): 

Proposition 7 Let g : R-|- R-|- be a monotone nondecreasing function. Con- 
sider a particle whose position changes at integer moments and is always an in- 
teger. If the particle is at position m > 1, it proceeds at the next .step to position 
m ~ X , where X is a random variable ranging over the integers 1, . . . , m — 1. 
All we know about X is that E[X] > g{m) and that X is chosen independently 
of the past. 

Let T be the random variable denoting the number of steps in which the 
particle reaches position 1. Then E\T\ < J" 

Applying this result to the potential function counting the number of unbal- 
anced triangles in graph T„, and employing Lemma [6] we infer that T5b(T„) = 
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6 The case p = 1/3 on a 2-regular triadic simplex 

Another case where the dynamics is easy to analyze is that when p — 1/3 and 
graph G is a 2-regular triadic simplex. In this case, according to Thcorem[TJ all 
recurrent states are either socially balanced (fixpoints), or all of them contain 
exactly one imbalanced triangle. Radicchi et al. [25] experimentally study the 
the convergence time on the triangular lattice and obtain the estimate 9{N°'), 
with a - 2.24. 

We can estimate the convergence time as follows: 

Theorem 4 For any 2-regular triadic simplex G such that Ty,{G) is connected 
there exists C > such that the time to social balance on G satisfies 

TuiTziG)) < Tsb{G) < TcRwiniG)) < C ■ „™™{3.2+log,(rc(T3(G))}^ (4) 

where tc{T3{G)) is the Cheeger constant of the triadic dual T3{G) and n is the 
number of vertices of T^{G) (i.e. the number 0/ triangles of graph G). 

Before going into the proof, note that an open problem in 2J (Chapter 14, 
Open problem 13), combined with Proposition 5 (same chapter) would imply 
that the lower and upper bounds in Theorem |4] stated in terms of graph T^{G) 
have the same order of magnitude. 
Proof. 

Since G is 2-regular, T^{G) is a 3-regular graphs has no loops, and the clique 
switching process reduces to an annihilating random walk. 

The annihilation time of an ARW is stochastically dominated by the coales- 
cence time of a CRW on the same graph. This is an easy example of a coupling 

m- 

Proposition 8 We can couple the annihilating and coalescing random walk 
such that Tarw Tcrw- Consequently 

Tarw{G) < TcRw{G). 

Proof. 

The coupling uses a well known idea (see e.g. [H]). We first describe the 
process intuitively. Consider a coalescing random walk. Declare particles that 
coalesce to be "ghost" particles (as opposed to the "live", not yet coalesced 
ones). If a set of ghost particle meet a real particle they stick to it, adopting the 
(randomly chosen) trajectory of the live particle. When two live particle meet 
they become ghosts. The coupling is obtained by seeing the ARW as the CRW 
restricted to live particles. 

Formally, define {S,t )x£V,t>o to be a stochastic process such that 

• Q = X for all X . 

• = implies that ^f+, = for all s > 0. 

• For all X G V, {£,t)t>o is a random walk on G. 
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Define the annihilating system rjt Q V, 

Vt^{y\ \{xeV,^^^y}\ is odd}. 

Assume that initially the number of particle was even. It is easy to see 
that when all particles have coalesced then all particles are ghost particles. The 
reason is simple: 

• Any cluster of particle contains at most one live particle. 

• The only way for a live particle to become a ghost is to meet another live 
particle. 

• The number of live particle only goes down by two at a time. 



For the upper bound we apply a result due to Aldous and Fill (Proposition 
9, Section 3.14, Chapter 14 in [2]), generalizing an earlier result of Donnelly and 
Welsh |lClj3 and infer the theorem. We use the fact that the lattice graph is 
3-edge-connected and 3- regular. 

For the lower bound we consider all initial states consisting of exactly two 
imbalanced triangles, and use the connection with annihilating random walks. 

❖ 

7 Application to 3-XOR SAT 

Our problem also displays an unexpected connection with a problem in the 
area of satisfiability solving: random 3-XOR SAT. The satisfiability of random 
instances of this problem has been investigated in both Statistical Physics [27] 
and Theoretical Computer Science [Hj. 

Definition 16 An instance F of the 3-XOR SAT problem is specified by a list 
of m equations on n variables, Xi-^ © Xi^ Xi^ = bi, for some bi £ {0, 1} and 
*i, J2, *3 G 1 . . . rt. F is satisfiable if there exists an assignment A of variables in 
F that makes every equation evaluate to true. 

Our results will not involve random instances, but satisfiable instances. 
Moreover, instances we will work with will be reduced, i.e. they satisfy the 
following two conditions: no variable appears in a single clause and no two 
clauses share two variables. This last assumption does not particularly con- 
strain the class of instances we want to solve. Indeed, it is easy to see that 
any formula can be transformed to a reduced one: We simply eliminate clauses 
involving pure variables (since they can always be satisfied by setting the pure 
variable the right way). Also, if Ci and C2 share two variables, say x and y, 
let z and t the remaining variables. The conjunction of Ci and C2 entails a 

^The result of Aldous and Fill is stated in [2 for continuous time, but can be easily 
translated to discrete time, at the expense of an additional linear factor 
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Algorithm RandomWalkSat($): 

Start with an arbitrary assignment U. 

while (there exists some unsatisfied clause) 
pick a random unsatisfied clause C 
change the value of a random variable of C in 

return assignment U. 



Figure 5: The Random WalkSat algorithm 

constraint of the type z ® t = A, for some A G {0, 1}. We can thus eliminate 
one of the variables z, t (by replacing it with t © A). This eliminates one of the 
clauses Ci and C2 from the formula. 

We will analyze the Random Walk algorithm displayed in Figure [S] We will 
assume that the input to the Random Walk algorithm is a reduced formula with 
the following additional properties: 

Definition 17 A 3-XOR formula is connected if one cannot partition F into 
two variable-disjoint formulas. F is fc-connected if one can delete up to k — 1 
variables (and the clauses of F where these variables appear) without discon- 
necting the formula. F is 2-regular if every variable appears exactly in two 
clauses. 

The displayed connection with Coalescing Random Walks enables us to prove 
the following result: 

Theorem 5 // reduced formula $ is satisfiable, has m equations, n variables, 
is 2-regular and s-connected then the expected time until RandomWalk finds a 
solution satisfies: 

E[Trw] <min{TOV2s,21og2-Tc(r3($)) -m^}. 

Proof. Consider the hypergraph s5(<i>) associated to the reduced formula and 
its triadic dual T<i,. Since $ is 2-regular it follows that T$ is actually a graph, 
and the dynamics of The RandomWalk algorithm on $ can be interpreted as 
an annihilating random walk on T$. We then apply the (discrete time version 
of) results of Donnelly and Welsh, as extended by Aldous and Fill in Chapter 
14 in [2]. 

8 Conclusions 

Our paper raises the interesting prospect that the convergence time of local 
search algorithms, such as the RandomWalk algorithm might be analyzed with 
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ideas from the interacting particle systems literature. We believe that this would 
be especially interesting for "message-passing" algorithms such as the belief and 
survey propagation algorithms [S] . 

Clearly, we would like to see progress in analyzing the dynamics of the ran- 
dom hyperedge switching process on hypergraphs. In particular studying the 
mixing time of the random hypergraph switching process is a very interesting 
problem for further study. Another possible direction, motivated by the analogy 
with the dynamics studied in |12j , is to further study the time to social balance 
using ideas similar to those in [23] . 

Finally, an interesting issue it to compute the Cheeger constant of the tri- 
angular lattice, as well as other planar regular lattices [15j . and comparing the 
results with those in [25] . 
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